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Water in the Near IR spectrum of Comet 8P/Tuttle 
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ABSTRACT 

High resolution spectra of Comet 8P/Tuttle were obtained in the frequency range 
3449.0-3462.2 cm^^ on 3 January 2008 UT using CGS4 with echeUe grating on 
UKIRT. In addition to observing solar pumped fluorescent (SPF) lines of H2O, the 
long integration time (152 minutes on target) enabled eight weaker H2O features to be 
assigned, most of which had not previously been identified in cometary spectra. These 
transitions, which are from higher energy upper states, are similar in character to the 
so-called 'SH' lines recorded in the post Deep Impact spectrum of Comet Tempel 1 
(Barber et al., 2007). We have identified certain characteristics that these lines have 
in common, and which in addition to helping to define this new class of cometary 
line, give some clues to the physical processes involved in their production. Finally, 
we derive an H2O rotational temperature of 62 ± 5 K and a water production rate of 
(1.4 ±0.3) X 10^8 molecules s'^ 
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1 INTRODUCTION 

Cometary nuclei are the most primitive objects in the solar 
system. Within about 3 AU of the sun, solar heating causes 
the surface temperature of the nucleus to rise above 200 K; 
water ice sublimes and other volatiles and icy dust grains 
are expelled from the surface. The motion, chemistry and 
rate of sublimation of molecules from the icy grains are in- 
fluenced by the solar wind and solar radiation. Initially the 
escaping gases form a dynamic, gravitationally-unbound at- 
mosphere, the coma. They are subsequently swept into the 
comet's plasma and dust tails and eventually dissipate into 
interplanetary space. The surface of the nucleus is processed 
during the course of repeated visits to the inner solar sys- 
tem, becoming covered with a rubble blanket of particles too 
large to be dragged away by escaping gas. During the course 
of numerous perihelia, the outer surface becomes completely 
depleted of volatile materials which either escape from the 
surface or else migrate inwards. This creates an outer man- 
tle of siliceous dust (Huebner, 2008), above a layered struc- 
ture in which the more volatile species are found at greater 
depths (Prialnik et al., 2008). This structure results in the 
great variety that is observed in active comets. 

Cometary nuclei typically have radii of ~3 km and are 
not able to be resolved from Earth. During close approaches, 
the inner coma of some comets are resolvable with large 
terrestrial instruments. However, in general, data obtained 
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with Earth-based instruments are indicative of conditions 
in a wider region of the coma, and our understanding of 
conditions within the inner coma (a region extending up to 
a few hundred km, where energy transfer is coUisionally- 
dominated, relies largely on modelling. Our understanding 
of the more extended regions of the coma, where the den- 
sities of the species are too low for thermodynamic equilib- 
rium to exist, also relies on modelling. In particular, mod- 
els of cometary coma predict the existence of so-called 'so- 
lar pumped fluorescent' (SPF) H2O emission lines. These 
originate from ro-vibrational excited upper states of the 
molecule, which, if observed in a coUisionally-dominated re- 
gion would be characteristic of kinetic temperatures of sev- 
eral thousand K, but are able to be produced in the cold, 
rarefied conditions existing in cometary coma, since these 
excited states, which have relatively large Einstein A coef- 
ficients, have time to decay radiatively before they are de- 
excited by collision with another molecule. 

Despite the success of models in predicting SPF lines 
(Crovisier, 1984; Weaver and Mumma, 1984, Mumma et al., 
1995), the physical conditions existing in cometary coma are 
highly complex and recent observations of Comet Tempel 1 
(Barber et al., 2007) revealed a type of emission line, not 
previously recorded in cometary spectra, and not predicted 
by the existing cometary models. These so-called 'SH' lines 
originate from upper states that are excited to higher vibra- 
tional energies than are the upper states of the SPF lines. 
Consequently, their production mechanisms are not able to 
be explained simply, by the existing cometary models. The 



2 R. J. Barber et al. 



near-Earth approach of Comet 8P/Tuttle in January 2008, 
provided an opportunity for us to search for, assign and char- 
acterise SH emission lines in its coma. We believe that our 
findings, which are presented in this paper, provide the basis 
for a more thorough investigation of the physical processes 
involved in the production of SH lines in cometary coma. 



2 COMET 8P/TUTTLE 

8P/Tuttle is a short period comet, Porb = 13.51 yr, and is 
the parent of the Ursid meteor stream (Jenniskens et al., 
2002). It had been estimated that the comet has a radius of 
~ 7.5 km, making it the largest of the group of 18 short pe- 
riod comets reported on by Licandro et al. (2000). However 
radar images obtained on 2-4 January 2008 (Harmon et al., 
2008) show a bifurcated nucleus consisting of two lobes each 
3-4 km in diameter, suggesting a possible contact binary. 
Photometric measurements, have yielded rotational times 
of the nucleus as being either 5.71±0.04 hr (Schleicher and 
Woodney, 2007b), 7.7±0.2 hr (Harmon et al., 2008), or mul- 
tiplicities of 5.7 and 7.4-7.6 hr (Drahus et al., 2008). How- 
ever, these results may need to be re-interpreted in the light 
of the probable bifurcation of the nucleus. 

In the scheme proposed by Levison (1996), comets are 
classified by reference to the Tisserand invariant. To a first 
order, this parameter is a constant of the motion of a comet. 
It is based on the Jacobi integral in the restricted three- 
body problem (Tisserand, 1894; Moulton, 1947). Because 
8P/Tuttle has a Tisserand invariant of less than 2 (Tj = 
1.601), it is classified as a 'near isotropic comet' (NIC) of the 
Halley type. However, its present orbit is not characteristic 
of this class, since most NIC comets have long orbital periods 
(defined as being greater than 200 yr). 

8P/Tuttle also differs from most short period comets as 
these normally have a Tisserand invariant 2 ^ Tj 5C 3 and in 
the Levison scheme are referred to as 'ecliptic comets' (EC). 

Being a short period NIC, the chemistry of 8P/Tuttle 
is of particular interest. Unlike most short period comets, 
which are of the EC type, 8P/Tuttle is believed to have been 
formed in the region of the giant planets. It therefore has the 
potential to reveal possible differences between EC and NIC 
objects. Because it has made many approaches to the Sun, 
8P/Tuttle's surface may be highly processed like that of an 
EC. If then, despite having a highly-processed surface, its 
spectra reveal a chemistry that is more like that of NICs 
than of ECs, it may be possible to identify differences in the 
chemistries of various regions of the very early solar system. 

Even though, with the exception of the 1953 approach, 
8P/Tuttle had been seen on each orbit since its discovery 
in 1858, until recently there had been no detailed investi- 
gation of its near-IR spectrum. During the 2008 apparition, 
8P/Tuttle approached to within 0.25 AU of the Earth and 
was favourably placed for observing (see for example, Bonev 
et al., 2008 and Bohnhardt et al., 2008). 

3 WATER IN COMETARY SPECTRA 

In the near IR, the spectra of comets are rich in ro- 
vibrational H2O emission lines which can, in principle, be 



observed from Earth. However, because of the low tempera- 
tures that characterise cometary comae, the strongest water 
lines are fundamental transitions (that is to say, transitions 
to ground vibrational states) from low-energy ro-vibrational 
states. Photons from these transitions are absorbed by wa- 
ter vapour in the Earth's atmosphere, the molecules of which 
are in ground vibrational states. 

Close to the nucleus the density is sufficiently great 
(generally upward of 10® molecules cm~"^) to ensure that the 
molecules are coUisionally thermalised in low J ground vibra- 
tional states, with a rotational temperature equal to the lo- 
cal kinetic temperature (Weaver and Mumma, 1984) . There 
is no precise boundary between the coUisionally-dominated 
and the less dense, radiatively dominated region, as the 
breakdown of thermal equilibrium occurs at different dis- 
tances from the nucleus for each of the allowed rotational 
transitions between ground vibrational states. These dis- 
tances vary between several tens of kilometres to several 
thousands of kilometres (Bockelee-Morvan, 1987). Moreover, 
beyond the region where H2O-H2O collisions are able, alone, 
to maintain local thermodynamic equilibrium (LTE) there is 
a region where the additional contribution of H2 0-electron 
collisions and line-trapping effects maintains LTE to greater 
radial distances (Bockelee-Morvan, 1987; Zakharov et al., 
2007). In particular, Xie and Mumma (1992) showed that 
at distances of several thousand kilometres, e~-H20 colli- 
sions play an important role in exciting water molecules. In 
SectionQwe discuss a number of mechanisms, including e~- 
H2O collisions, that may possibly be involved in the forma- 
tion of the so-called 'SH' water lines that were observed to 
be present in our 3 January 2008 UT spectrum of 8P/Tuttle. 
SH, or 'stochastic heating' lines were first identified in the 
'Deep Impact' spectrum of Comet Tempel 1 (Barber et al., 
2007), and were subsequently noted to be present in another 
spectrum of the same comet (Mumma et al., 2005, Figure 2). 
Their characteristics are discussed in Section[5] 

At distances greater than a few thousand kilometres 
from the nucleus the density of the species is sufficiently low, 
and the mean free path between collisions large enough, for 
'fluorescence equilibrium' to exist. That is to say, radiation is 
due to the balance between solar pumping to vibrationally- 
excited states and subsequent spontaneous radiative de- 
cay (Crovisier, 1984). Under coUisionally-thermalised con- 
ditions, these upper states would typically require tempera- 
tures of several thousand Kelvin for them to be sufficiently 
populated for their coUisionally-induced radiation to lower 
states to be observable. The observation in cometary spectra 
of transitions from vibrationally-excited upper states, indi- 
cates that they must be produced in regions where Boltz- 
mann statistics do not apply. 

When vibrationally excited states decay radiatively to 
lower energy ro-vibrational states above the ground vibra- 
tional state (hot bands), the photons emitted are not ab- 
sorbed in the Earth's atmosphere. Recognition of this fact 
was a key factor in developing the 'solar-pumped fluorescent' 
(SPF) approach to analysing cometary spectra (Mumma et 
al., 1995). Spectra containing these SPF transitions enable 
information to be derived about the rotational temperature. 
Trot of the inner coma and the water production rate, Qh2 0- 
A comparison of this latter value with the production rates 
of trace species provides the abundance ratios that are fun- 
damental to understanding cometary chemistry. 
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Because the upper states populated by solar pumping 
are not coUisionally thermalised, the derivation of physical 
quantities from line intensities requires knowledge of all the 
possible excitation as well as the alternative de-excitation 
routes. Each upper excited ro- vibrational state can be popu- 
lated directly by radiative pumping of rotational levels of the 
ground vibrational state, as well as by cascade from higher 
energy states (Crovisier, 1984). The relative probabilities of 
these various upward and downward routes are temperature- 
dependent. Models have been developed to compute the g- 
factors for each SPF line, from which line intensities are able 
to be calculated as a function of temperature (Crovisier, 
1984; Mumma et al., 1995, Dello Russo et al., 2000). These 
models require a knowledge of the Einstein A coefficients 
for the downward transitions. Ideally, the upward pump- 
ing rates should be computed using the Einstein B coeffi- 
cients. However, for simplicity, the temperature-dependent 
vibrational band intensities are frequently used to compute 
the upward pumping rates. For example, Dello Russo et al. 
(2004) used Einstein A coefficients from the BT2 synthetic 
water line list (Barber et al., 2006) to model 5- factors for all 
the H2O SPF transitions in the: vi-\-Vi-i'\, I'l+vs-i'a, 2u\- 
u\, 2vi-vz and v\-\-i/2-\-i"j,-i'i-V2 vibrational bands, from 
upper states having J ^ 7. 

The 2.9-3.0 fim (3300-3450 cm"^) spectral region is 
particularly rich in SPF lines and is also largely devoid 
of other molecular species. The availability of accurate g- 
factors has enabled this region to be used in determin- 
ing cometary rotational temperatures and water production 
rates (e.g. Dello Russo et al., 2004 and 2006). 



remove the sky. Each ABBA series represented a total of 
160 seconds integration time. Total time on-target was 152 
minutes. Spectra of an AOV star, HD 6457, were obtained 
at the beginning and end of the observing session. Standard 
star spectra are required in order to adjust for frequency- 
dependent atmospheric transmission, and also to enable ab- 
solute flux calibration. Before dividing by the standard star, 
we removed high frequency noise from both the 8P/Tuttle 
and standard star data using Fourier transform smoothing. 
We rejected all data at frequencies where the atmospheric 
transmission was less than 20%, and this resulted in 19% of 
the data in the 3 440.6 to 3 462.6 cm~^ range being excluded. 
We noted that transmission was poor at frequencies less than 
3449.0 cm~^ and above 3 460.2 cm~^. We therefore decided 
to limit the spectral region under investigation to that ly- 
ing between these frequencies. In this reduced, 11.2 cm~^, 
frequency range less than 3% of the data failed to meet our 
20% minimum atmospheric transmission requirement, and 
the average transmission in this reduced frequency range 
was 66%. After dividing by the standard star, flux calibra- 
tion was achieved by multiplying the result by the frequency- 
dependent flux of HD 6457, calculated relative to Vega as 
a photometric 0.0 magnitude standard (Hayes, 1985), and 
assuming a black body function with Tcff — 9 480 K (Toku- 
naga, 2000, pl51). 

Because of the lack of arc lines in this spectral region, we 
frequency-calibrated our data using the position of several 
known SPF lines that were present in the 8P/Tuttle spec- 
trum. This had the effect of transforming the observed spec- 
trum into the rest frame. The fully-reduced, flux-calibrated 
spectrum is given in Figure 1. 



4 OBSERVATIONS 

We observed 8P/Tuttle on its approach to the Sun (peri- 
helion being on 27 January UT), on the nights of 3, 4 and 

5 January 2008 UT from UKIRT, Hawaii, using the echelle 
grating on CGS4. We report here on the spectra that we 
obtained on 3 January UT, in the frequency range 3 440.6- 
3 462.6 cm~^. The heliocentric distance on this date was 
~1.09 AU, and the geocentric distance and velocity were 
~0.25 AU and ~3.3 km s~^ respectively; the latter fig- 
ure corresponds to a red shift of 0.038 cm~^, which was 
equivalent to 1 pixel on the array, and less than the mini- 
mum resolution of our instrument, which at R = 37 000 was 
0.093 cm~^ at 3 450 cm~^. The total visual magnitude of 
8P/Tuttle at the time of our observations was ~5.8. How- 
ever our observations were of the bright inner coma region 
which presented a point source in our instrument, and re- 
quired long integration times. 

We used a slit of 2-pixel width of length 90 arcsec, ori- 
ented east-west. The scale was 0.86 arcsec per pixel in the 
spatial direction (156 km at the comet) and 0.41 arcsec in 
the spectral direction, (149 km at the comet for a 2-pixel 
slit). The telescope was dithered to give 2x2 sampling ev- 
ery 40 seconds. Observations were acquired using a standard 
ABBA sequence, which was achieved by nodding the tele- 
scope 22 pixels along the slit. In this mode, the comet signal 
was present in both the A and B beams (rows 114 and 92 
respectively), which, compared to nodding to blank sky, in- 
creased the signal to noise ratio by a factor of \/2, when 
both signals were added using the (A-B)-(B-A) routine to 



5 DATA ANALYSIS 

We compared the features in Figure 1 with our BT2 
database of H2O SPF line positions and factors. Nine of 
these features are at frequencies corresponding to known 
SPF lines (there are actually eleven SPF lines, with two fea- 
tures being blends of two SPF lines). Details are given in 
TableO However, Figure 1 also includes a number of other 
features, whose shape and intensity above the continuum 
(S/N ^2) suggest that they are spectral lines, rather than 
noise, but whose frequencies do not correspond to SPF lines. 

In order to assign these other features, we employed the 
methodology detailed in Barber et al. (2007) for the identi- 
fication of non-SPF features in the post Deep Impact spec- 
trum of Comet Tempel 1. Initially we used the BT2 line list 
to generate synthetic emission spectra for water at a tem- 
perature of 3 000 K in the 3 449.0-3 462.2 cm"^ frequency 
range, restricting the value of J for the upper state to 8, and 
the intensity of the weakest line to be generated was lim- 
ited to 10~^ of the strongest. The output, which consisted 
of about 600 individual lines was convolved to match the 
resolving power of the CGS4 echelle, R=37000. 

Some of the features in this synthetic spectrum are at 
frequencies corresponding to features in the spectrum of 
8P/Tuttle. However, the relative intensities of the features 
in the observed and synthetic 3 000 K spectra do not agree. 
The reason is that the observed emission lines originate in 
low-density regions where vibrationally-excited upper states 
have sufficient time to decay radiatively before they can be 
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Figure 1. The upper chart shows the spectrum of Comet 8P/Tuttle, UT 3 January 2008 adjusted to rest frame and plotted after dividing 
by standard star HD 6457 (normalised) to adjust for frequency-dependent atmospheric transmission. The spectra of both the comet and 
standard star have been Fourier Transform smoothed to remove high frequency noise. Atmospheric transmission is shown in the lower 
chart, and no spectral data is displayed where transmission is less than 0.20. The figures above the spectrum indicate lines that have 
been assigned in Table[2] the numbering above the plot corresponds to that given in each entry in this table. Bold font is used to indicate 
SH transitions; SPF lines are numbered using normal weight font. The letters below the plot correspond to the nine SH lines which we 
would expect to be present in the spectrum (see text), but which were not positively identified. Their numbers correspond with those 
given in the first column of each entry in Table[3] Six of these either correspond to weak feature (S/N 2), or are at frequencies where 
they would be undetectable due to blending with stronger lines. Of the twenty-six SH lines listed in TableO we are unable to explain 
the absence of three, and these are identified with bold letters. 



coUisionally de-excited. These upper states are populated 
by pumping low-lying rotational states and by cascade from 
higher levels (of which there will be many, most of which 
contribute little and can therefore be disregarded). Since the 
water molecules are not coUisionally thermalised, the popu- 
lations of the upper states can be substantially greater than 
would be the case if Boltzmann statistics applied, and the 
frequencies and intensities of the emissions from these states 
are frequently similar to those of emission lines from H2O 
vapour in LTE regions at temperatures of, say, 3 000 K. 

Our synthetic spectrum contains many more features 
than are present in the spectrum of 8P/Tuttle. It is there- 
fore important to try to identify specific characteristics of 
the observed lines that differentiate them from lines in the 
synthetic spectrum that are absent in the observed spec- 
trum. 

When synthetic spectra are generated using BT2, one 
of the output files gives the frequency, assignment, inten- 
sity, Einstein A coefficient and lower state energy of every 
water transition in the frequency range, down to a specified 
intensity cut-off. Using these data, we were able to identify 



certain characteristics of those lines in our synthetic spectra 
whose positions matched those of features in the observed 
spectrum, that differentiated them from those that did not. 
This information enabled us to constrain the parameters ap- 
plied when generating the BT2 synthetic spectrum, which 
reduced the number of lines in the synthetic spectrum and 
hence the risk of confusing noise and real spectral lines. 

Before detailing these tighter parameters and our find- 
ings, we comment briefly on the ro- vibrational states of wa- 
ter. The three vibrational modes of the H2O molecule, vi, 1/2 
and !/3 (symmetric stretch, bend and asymmetric stretch 
respectively), require different amounts of energy for one 
quanta of excitation. Also, the magnitude of each increment 
decreases as the total internal energy of the excited molecule 
increases. In the energy range in which we are interested, 
each quanta of symmetric and asymmetric stretch represents 
~3 400-3 700 cm~^. These values are approximately twice 
the energy of one quanta of bend. Hence it is convenient to 
talk in terms of 'polyads', where each polyad, designated u 
is the equivalent of one quanta of stretch or two of bend. 
Single bend quanta are denoted by '5'. 
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We noted that all the SPF lines in the observed spec- 
trum came from the 2v polyad, whilst the large majority of 
non-SPF lines that matched features in our synthetic spec- 
tra, were from the 'iv and <iv+5 polyads. We also noted 
that none of the upper states had more than two V2 quanta. 
In addition, the lines in our synthetic spectrum that corre- 
sponded to non-SPF features in Figure 1 all had Einstein A 
coefficients J5 f .0 s~^, and were on average higher than those 
of the SPF lines, the average A^/ being 23.0 and 7.4 s~^ re- 
spectively. The SPF lines had upper states with energies in 
a very narrow range: 7 242-7 613 cm~^, whilst the non-SPF 
lines had upper states with energies in the range 10 365 - 
12 940 cm-\ 

Rotational excitation, also affects which lines are ob- 
served in cometary spectra. This is defined in terms of the 
asymmetric top labels, J, K^, Kc which represent total an- 
gular momentum and the prolate and oblate levels respec- 
tively. We designate ro- vibrational states {vi V2V3)[J K^, ifj 
and when assigning transitions, give the upper state first. 

An examination of the ~ 600 transitions in our syn- 
thetic spectrum revealed that none of the features in the 
spectrum of 8P/Tuttle corresponded to lines whose upper 
state had J ^ 5. We also noted that none of the lines in 
our observed spectrum, assigned in Tabled correspond to 
very weak features in our synthetic spectrum. They all cor- 
respond to lines in our synthetic spectrum having intensities 
greater than ~0.5% of the intensity of the strongest line in 
the synthetic spectrum. In all of these respects the non-SPF 
lines that we identified in the spectrum of 8P/Tuttle had the 
same characteristics as the SH lines observed in the 'Deep 
Impact' spectrum of Comet Tempel 1 (Barber et al., 2007). 

This information provided us with with a second, even 
tighter, set of parameters with which to search for SH lines. 
We produced a second synthetic BT2 spectrum with the 
constraints: J ^ 2>v or 2>v-'rS polyad, ^if ^ 1.0 s~^ and a 
minimum intensity cut-off equivalent to 0.5% of the intensity 
of the strongest line in the frequency range. Also we lowered 
the temperature of the synthetic spectrum to 2 500 K, which 
further reduced the total amount of data generated. 

Although it still includes all the SPF lines that were 
present in our observed spectrum, this second synthetic LTE 
spectrum contains only 26 possible SH lines. Using these 
data, we were able to assign seven unblended SH lines in 
the observed spectrum. These are detailed in Tabled which 
also gives Euppor, A^/; the type of transition, SPF or SH; 
the nuclear spin species identity (ortho/para) and the esti- 
mated signal-to-noise ratio of the lines. In addition to the un- 
blended SH lines, we identified a feature centred at 3458.49 
cm-^ There is a known SPF line, (200)[423]-(100)[432], at 
3458.51 cm~^. However, at 62 K, the p-factor for this line is 
only 8.2x10"^", which equates to a line of only ~3% of the 
intensity of the feature observed at 3 458.49 cm~^. The BT2 
synthetic spectrum contains only one other water line in the 
3458.49 cm"^ region, an SH line, (022)[432]-(120)[541] at 
3458.47 cm~^. We therefore conclude that this line is also 
present in the observed spectrum of Comet 8P/Tuttle, which 
represents an eighth SH detection. 

Barber et al. (2007) assigned two unblended SH transi- 
tions in the spectrum of Comet Tempel 1 at 3 453.39 and 
3 451.51 cm~^ and a feature at 3453.90 cm~^ which they 
suggest might be a blend of two lines: (211)[322]-(210)[221], 
which has a ?iv~\~5 polyad upper state, of energy 12 354 cm~^ 



and (103)[110]-(102)[110], which has a 4,v polyad upper state 
of energy 14356.8 cm"^ The spectrum of 8P/Tuttle in Fig- 
ure 1 also contains a feature at 3 453.90 cm~^, and based on 
the polyad constraints that we detail above, we are inclined 
to believe that this is the unblended (211) [322]-(210) [221] 
line, and that there is no observable (103)[110]-(102)[110] 
component. This view is reinforced by considerations of K^. 
and (see below). Figure 1 also contains a feature (S/N = 
2.1) at 3 451.51 cm"\ which we identify as the (220) [212]- 
(021)[111] SH transition, which is in agreement with Barber 
et al's assignment. However, there is no feature in our ob- 
served 8P/Tuttle spectrum at 3453.39 cm~^. This may be 
because of the width of the strong feature (a blend of two 
SPF lines) centred at about 3 453.25 cm~^, since a line at 
3 453.39 cm~^ would be in wing of this feature. It should 
also be noted that the line that Barber et al identified in 
Tempel 1 at 3453.39 cm"^ is the (210)[101]-(011)[110] tran- 
sition, whose upper state is of the 2^+6 polyad, at an en- 
ergy of 8 784.7 cm"'^, and hence would not satisfy the 3u and 
3u+S test that we are applying for SH lines in the current 
paper. 

Having identified that the observed spectrum of 
8P/Tuttle contains eight of the 26 SH fines appearing in 
TableO we examined the data for additional characteristics 
that might distinguish the observed eight from the unob- 
served 18, and noted one further feature that is common 
to all the observed SH lines; they all have Ka, J5 J-1. Now 
of the 18 lines in Table 3 that were not observed, half had 
Ka, ^ J-1, and half did not. Moreover, of the nine lines 
in Tab leO that had A'a ^ J-1, but which were not identi- 
fied in the observed spectrum, only the (300)[542]-(101)[441] 
transition at 3 457.85 cm~^, labelled b in Table|3]was defi- 
nitely absent. Adopting a cautious approach we also include 
the (220)[322]-(120)[413] transition at 3 452.42 cm"^ and the 
(003)[440]-(002)[541] transition at 3 449.64 labelled m and 
r respectively in TableO as being not present even though 
there is a weak feature in Figure 1 centred at 3 452.48 cm~^, 
which is only 0.06 cm~^ away from m, and therefore theo- 
retically unresolvable from it, and a 2.1 S/N feature centered 
at 3 449.55, which would also theoretically be unresolvable 
from the synthetic SH line r. Of the remaining six lines in 
TableOthat have Kn, ^ J-1, but were not observed to be 
present in Figure 1, two are at frequencies where they would 
be blended with stronger SPF features, and four are at fre- 
quencies corresponding to weak features (S/N 5C 2). On this 
basis, we estimate the probability as 1 in 763 that the prop- 
erty common to all SH lines assigned here (JCa ^ J-1) is a 
chance event. We therefore conclude that this characteristic 
has a physical cause, most probably relating to the mech- 
anism by which SH lines are produced. We discuss this in 
more detail in Section[71 but comment here that for a given 
J, there are 2J+1 possible combinations of and Kc, and 
the higher the value of Ka_, the higher is the energy of the 
state. The preferential population of higher value JsTa states 
is the reverse of what would be observed in a Boltzmann dis- 
tribution, and this inversion of states strongly suggests that 
the upper levels of SH transitions are populated by cascade 
from more energetic levels, rather than by pumping from 
ground vibrational states, which is the principal method of 
populating the upper states of SPF transitions. 

In addition to analysing our own 8P/Tuttle data, we ex- 
amined the spectra of comet C/1999 HI Lee obtained on 19 
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and 21 August 1999 using NIRSPEC on Keck (Dello Russo 
et al., 2006). These contained several weak, unassigned, fea- 
tures. We noted that the frequencies of three of these fea- 
tures corresponded to SH lines in the spectrum of 8P/Tuttle. 
These are: 

(022)[431]-(120)[542] at 3459.85 cm"^ 
(211)[322]-(210)[211] at 3 453.90 cm"^ 
(220)[212]-(021)[111] at 3451.51 cni'V 

We do not claim that our assignments are totally secure. In 
particular we observe that the feature at 3 451.51 cm^^ was 
only present in the spectrum obtained on 19 August, whilst 
a line at 3 453.90 cm~^ would be likely to appear blended 
(see below). However, the fact that these two lines were 
also observed in Comet Tempel 1 increases our confidence in 
these assignments. The third line, at 3 459.85 cm~^ is out- 
side the frequency range examined in Barber et al. (2007). 
Dello Russo et al. state that the S /N ratio of the feature cen- 
tred at ~3 453.85 cm~^ in their spectra of comet Lee is 7.9, 
and indicate that its position corresponds to the (110) [313]- 
(010) [422] SPF line at 3 453.88 cm~\ However, the g factor 
for this line is too low to account for a signal of this intensity, 
which leads us to believe that the main component in the 
blend is likely to be the (211)[322]-(210)[211] SH line. Com- 
paring the intensities of the features at 3 459.85 and 3 451.51 
cm~^ on the two nights with that at 3453.85 cm~^, we in- 
fer that the signal-to-noise ratios of the SH lines alone are 
between 2.5 and 3.5. 

Using the original data file and applying Fourier Trans- 
form smoothing to remove high frequency noise, we identi- 
fied several water lines that were unassigned in the high 
resolution spectrum of 8P/Tuttle obtained at ESO with 
CRIRES on 27 January 2008 UT, when the heliocentric dis- 
tance was 1.03 AU (Bohnhardt et al., 2008, Fig. IC). These 
include SPF lines at 3 451.09, 3 456.45 and 3 458.12 cm-\ aU 
of which are also present in Tabled We also noted that the 
spectrum contained an unassigned feature at 3 455.48 cm~^ 
which is the frequency of the (013) [330]-(012) [413] SH line 
in Table [2l We estimate the S /N ratio of this feature to be 
~3 and note that its width suggests that it may be a blend 
of more than one line. 

Finally, we note that the unassigned feature at 
3 453.90 cm~^ in the spectrum of Comet 73P/Schwassmann- 
Wachmann 3C obtained by N. Dello Russo on 15 May 2006 
UT (private communication), corresponds to the (211)[322]- 
(210) [211] SH line in Tabled 

On the basis of our identification of SH lines in comets 
of different types, we conclude that this class of water lines, 
which until recently had not been identified in cometary 
spectra, is probably always present at heliocentric distances 
of up to ~1.5 AU. However, it is noted that the SH lines 
observed vary from comet to comet, and indeed between 
spectra of the same comet obtained on different dates. Ta- 
ble[T] contains details of the four comets in which SH lines 
have been recorded. 



6 ROTATIONAL TEMPERATURE AND H2O 
PRODUCTION RATE 

The procedure that we used to determine the temperature of 
the inner coma is based on the fact that the p-factors of some 
SPF transitions can behave quite differently as a function 



Table 1. List of comets in which SH water lines have been ob- 
served 



Name 


P (yr) 


Type 


Origin 


rh 


T,ot K 


9P/Tempel 1 


5.5 


lEC 


Kuiper 


1.51 


~40 


C/1999 Hl(Lee) 


78 000 


NIC 


Oort 


1.06 


~78 


73P/SW3 


5.4 


2EC 


Kuiper 


1.03 


~110 


8P/Tuttle 


13.5 


3nic 


Oort 


1.09 


~62 



EC; Ecliptic Comet 

NIC: Near Isotropic Comet 

^ artificial impact 

^ fragmented 

^ bifurcated nucleus 



of temperature: some (^factors increase with temperature, 
whilst others decrease (Dello Russo et al., 2004). In broad 
terms, SPF lines whose upper states have been pumped from 
higher J values become more intense as temperature rises, 
whilst those pumped from the lower J states show little in- 
crease in intensity at higher temperatures, and in the case 
of those pumped from the lowest rotational states, may ac- 
tually weaken with increasing temperature. 

Using the g-factors for the SPF lines in our spectral 
frequency range (Dello Russo et al., 2004), we gener- 
ated SPF synthetic spectra, convolved to the resolving 
power of our instrument. The spectra were produced 
using the Fortran program, spectra-BT2 (Barber et 
al., 2006) which is available in electronic form via: 
http:/ / www. tampa.phys. ucl. ac. uk/ftp / astrodata /water / B T2\ 
and were generated at various temperatures with 5 K 
increments (the smallest temperature difference that pro- 
duced measurable differences in the synthetic spectra). The 
spectral region contains many SPF transitions in the five 
vibrational bands, for which ^-factors were computed by 
Dello Russo. Some of these lines are not detectable at low 
rotational temperatures, whilst others are at frequencies 
that fall within regions of low atmospheric transmission. 
We assigned nine SPF lines in Figure 1, of which two 
are blended. With the exception of the SPF/SH blend at 
3458.49 cm~^ (which as discussed above is predominantly 
an SH line) we used all of these features for temperature 
diagnostics. 

On comparing the relative intensities of the synthetic 
spectra with our observed spectrum, normalised to the 
cometary continuum level, we observed that there was excel- 
lent agreement between the relative intensities of all the SPF 
lines at 62 ±5 K, and consequently we conclude that this was 
the rotational temperature of the inner coma of 8P/Tuttle 
at UT 3.30 January 2008 (the mid-point of our observations) 
when the geocentric distance was ~0.25 AU. In arriving at 
this temperature, we have assumed the normal ortho-para 
ratio (OPR), which is 3:1. Some comets have been observed 
to have sub-normal OPRs. However, as all seven features 
used to derive the temperature of the inner coma are tran- 
sitions between ortho states, our results are not sensitive to 
the OPR. Our derived rotational temperature is consistent 
with those obtained by Bonev et al. on 22 and 23 Decem- 
ber 2007, of 60±15 K and 50±10 K respectively, when the 
geocentric distances were ~0.32 and ~0.31 AU respectively. 

Due to slit losses, nucleus-centred spectra provide a wa- 
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Table 2. Assignments of SPF and SH lines in the post-impact spectrum of 8P/Tuttle. In order, the columns give: a reference number: 
the rest frequency in cm"^; the assignment, upper level given first (vibrational quantum numbers in round brackets, rotational quantum 
numbers, J, K^, Kc in square brackets); energy of the upper level in cm~^; Einstein A coefficient from BT2; type of line, SPF/SH; 
ortho/para state; S/N ratio and comments 



Ref. 


Preq. 


Identification 




A,f 


Type 


O/P 


S/N 


Comment 




cm~^ 


see text 


cm ^ 


s-i 










1 


3459.85 


[ VZZ ) [4oiJ- 




10 934 


3.2 


SH 


P 


3.0 


Also seen in Lee^ 


2 


3459.53 


(101)[111]- 


(001)[202] 


7 285 


1.1 


SPF 


O 


4.2 


Possibly blended with 3, and a in Table 3 


3 


3459.49 


(101)[431]- 


(100) [532] 


7 613 


22.2 


SPF 


o 


4.2 


Possibly blended with 2 and a in Table 3 


4 


3458.51 


(200) [423]- 


(100) [432] 


7489 


1.1 


SPF 


o 


5.1 


Intensity suggests blending with 5 


5 


3458.47 


(022) [432]- 


(120)[541] 


10 933 


3.1 


SH 


o 


6.1 


Blended with 4, see above comment 


6 


3458.12 


(101)[000]- 


(001)[111] 


7 250 


3.6 


SPF 


o 


4.7 




7 


3458.00 


(300) [541]- 


(101) [440] 


11166 


4.9 


SH 


o 


4.7 




8 


3457.53 


(003) [432]- 


(002) [533] 


11382 


71.2 


SH 


p 


4.3 




9 


3456.45 


(101) [422]- 


(100) [523] 


7552 


31.9 


SPF 


o 


2.7 




10 


3455.98 


(200) [303]- 


(100) [414] 


7334 


4.5 


SPF 


o 


5.5 


Possibly blended with / in Table 3 


11 


3455.48 


(013) [330]- 


(012) [431] 


12 940 


51.3 


SH 


p 


3.5 


Also seen in 8P/Tuttlet 


12 


3454.69 


(101)[211]- 


(001) [220] 


7342 


1.8 


SPF 


o 


5.9 




13 


3453.90 


(211)[322]- 


(210)[221] 


12 354 


8.5 


SH 


o 


3.0 


Also seen in Leet, 8P/Tuttle*, Tempel 1 * 


14 


3453.30 


(200)[110]- 


(100) [221] 


7 242 


4.7 


SPF 


o 


13.3 


Possibly blended with i in Table 3 


15 


3453.15 


(101)[202]- 


(100)[321] 


7318 


1.7 


SPF 


o 


7.4 




16 


3451.80 


(121)[322]- 


(120) [423] 


10 550 


37.2 


SH 


o 


2.1 




17 


3451.51 


(220) [212]- 


(021)[111] 


10 365 


4.4 


SH 


o 


2.1 


Also seen in Leet SP/Tuttle*® Tempel 1* 


18 


3451.09 


(101) [413]- 


;001)[422] 


7517 


1.5 


SPF 


o 


2.3 




19 


3450.29 


(200) [110]- 


(001) [111] 


7242 


6.6 


SPF 


o 


14.6 


Possibly blended with o,p in Table 3 



t Dello Russo et al., 2006 * Bohnhardt et al., 2008 * Barber et al, 2007 ® marginal detection only 



ter production rate that is less than the global value derived 

from line intensity measurements beyond the seeing disk 
(Dello Russo et al., 2000). Hence, our derivation of the global 
production rate of gaseous H2O molecules in 8P/Tuttle in- 
volved an adjustment to compensate for slit losses. These 
were estimated by reference to the measured percentage of 
the signal diffracted into rows on the array (spatial direc- 
tion) that were adjacent to the row on which the image 
was focussed. We adopted a model in which all H2O is pro- 
duced close to the nucleus and released symmetrically into 
the coma with uniform velocity (Dello Russo et al., 2004). 
The result, based on the 3 Jan 2008 UT observations of 
8P/Tuttle, was 1.4 ± 0.3 x 10^* molecules s"\ 

There are several other estimates of the H2O produc- 
tion rates. Two are based on OH production rates ob- 
tained using narrow-band photometry. These are: 0.39 x 10^* 
molecules s^^ on 1 Nov 2007 UT, at a heliocentric distance 
of 1.63 AU (Schleicher, 2007) and 0.76 x 10^*^ molecules s"^ 
heliocentric distance of 1.30 AU, averaged from observations 
on 3, 4, 5 Dec 2007 UT, (Schleicher and Woodney, 2007a). 
Neither of these figures comes with an estimate of error. 
Bonev et al., (2008) give two values obtained on 2007 De- 
cember 22-23 UT, using different Keck NIRSPEC settings, 
when the comet was at a heliocentric distance of ~1.15 AU. 
These are: (2.38 ± 0.06) x 10^* and (2.13 ± 0.11) x 10^* 
molecules s^^. Bohnhardt et al., 2008, give water produc- 
tion rate of (5.97 ± 0.27) x 10^* on 27 January 2008 UT, 
when the heliocentric distance was ~1.03 AU. 

An examination of these estimates suggests that after 
taking account of the differences in the heliocentric distance 
of 8P/Tuttle on the various dates, our estimated production 
rate of 1.4±0.3 x 10^* molecules s~^ is low in comparison to 
other estimates, but is not inconsistent with them as H2O 



production rates in comets can vary considerably over short 
periods of time. 



7 ORIGIN OF SH LINES 

We comment here on the possible origins of the SH lines 
that we have identified in cometary spectra and which are 
characterised by upper states of the Sv and 3u+S polyad 
with excited symmetric/asymmetric stretch modes coupled 
with low bending (1^2 ^ 2) and and low rotational (J ^ 5) 
excitation, with only transitions from the more energetic 
Kg, > J — 1 states being observed. 

Population inversion {Ka, > J — ^) suggests that the 
upper states of the SH transitions are populated by cas- 
cade from higher vibrationally-excited states, rather than 
by pumping from lower levels (which would preferentially 
populate high K^, low if a states which are at lower en- 
ergies (see for example, the three para states of water: 
(300) [440], (300) [422] and (300) [404] which have energies, 
11048.4, 10898.1 and 10810.2 cm"\ respectively). The lad- 
der of downward transitions from an over-populated upper 
state through progressively lower energy levels is a 'water 
maser backbone series', and is illustrated in Figure 1 of 
Cooke and Elitzur (1985). The manner in which the high- 
est level of the backbone is populated is an area for fur- 
ther research. Here we mention three possibilities: (i) di- 
rect excitation by electrons (ii) dissociative recombination 
of HaO"*", and (iii) photo-dissociation of water with subse- 
quent quenching of 0(^D) by an H2O molecule. 

Xie and Mumma (1992) demonstrate that in the case 
of active comets, electron-water collisions can play an im- 
portant role in populating the rotational states. They base 
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Table 3. List of all the SH lines in a BT2 synthetic emission spectrum at 2500 K in the frequency range 3449.0 to 3460.2 cm"^, 
satisfying the empirical sort parameters detailed in the text. The standard convention is used for assigning the transitions (see Table 2). 
Lines that have been observed carry the same reference number as in Table 2, whilst those that arc not observed in Figure 1 are denoted 
by a lower case letter in the first column. It is noted that no observed line had J—Kai^ 1. Applying this additional sort parameter, we 
conclude that there arc a maximum of three SH lines which we might expect to be observed in this frequency region that are able to be 
identified as not being present in Figure 1. Units of intensity are Wm^-^cm mol~^. 



Ref. 


Freq. 


Identification 




O/P 




l2500/f 


J~Ka^ 1 


Comment 




cm~-^ 




see text 


cm ^ 




s-1 








1 


3459.85 


(022) 


[431]- 


(120)[542] 


10 934 


P 


3.2 


3.1E-19 


Yes 


Seen 


a 


3459.50 


(201) 


[514]- 


(101) [423] 


10 996 


P 


5.9 


6.7E-19 




Would be blended with 2 in Table 2 


5 


3458.47 


(022) 


[432]- 


(120)[541] 


10 933 


O 


3.1 


8.9E-19 


Yes 


Seen 


7 


3458.00 


(300) 


[541]- 


(101) [440] 


11 166 


O 


4.9 


1.5E-18 


Yes 


Seen 


b 


3457.85 


(300) 


[542]- 


(101) [441] 


11 165 


p 


4.9 


5.0E-19 


Yes 


Not seen and J—Ka^ 1 


8 


3457.53 


(003) 


[432]- 


^^^^^^^ rc ool 

002) [533] 


11382 


p 


71.2 


5.3E-18 


Yes 


Seen 


c 


3457.31 


(211) 


[321]- 


(210) [220] 


12 360 


p 


14.3 


4.7E-19 


Yes 


Coincides with a weak feature in Fig. 1 


d 


3457.28 


(220) 


[202]- 


(021)[101] 


10 352 


p 


6.6 


4.9E-19 






c 


3456.53 


(013) 


[321]- 


(012) [422] 


12 772 


p 


93.6 


2.4E-18 


Yes 


Coincides with a weak feature in Fig. 1 (blend) 


f 


3455.96 


(201) 


[202]- 


(200)[101] 


10 681 


o 


20.4 


3.8E-18 




Would be blended with 10 in Fig. 1 


11 


3455.48 


(013) 


[330]- 


(012)[431] 


12 840 


p 


51.3 


1.3E-18 


Yes 


Seen 


13 


3453.90 


(211) 


[322]- 


(210)[221] 


12 354 


o 


8.5 


8.4E-19 


Yes 


Seen 


g 


3453.51 


(310) 


[523]- 


(111)[422] 


12 575 


o 


15.0 


2.0E-18 






h 


3453.45 


(003) 


[422]- 


(002) [523] 


11332 


o 


84.7 


1.9E-17 






i 


3453.33 


(003) 


[431]- 


(002) [532] 


11384 


o 


73.2 


1.6E-17 


Yes 


Would be blended with 14 and 15 in Fig. 1 


j 


3452.90 


(003) 


[202]- 


(002) [321] 


11100 


o 


4.7 


6.8E-19 






k 


3452.63 


(21i; 


[303]- 


(210) [202] 


12 282 


p 


24.4 


8.4E-19 






1 


3452.42 


(300) 


[524]- 


(101) [423] 


10 989 


p 


6.9 


7.9E-19 






m 


3452.42 


(220) 


[322]- 


(120) [413] 


10 510 


p 


17.0 


1.6E-18 


Yes 


Not seen and J—Ka^ 1 


16 


3451.80 


(121) 


[322]- 


(120) [423] 


10 550 


o 


37.2 


l.OE-17 


Yes 


Seen 


17 


3451.51 


(220) 


[212]- 


(021)[111] 


10 365 


o 


4.4 


9.7E-19 


Yes 


Seen 


n 


3450.82 


(201) 


[212]- 


(200)[111] 


10 688 


p 


18.0 


l.lE-18 


Yes 


Coincides with a weak feature in Fig. 1 


o 


3450.24 


(211) 


[322]- 


(111)[211] 


12 354 


o 


6.0 


6.0E-19 


Yes 


Would be blended with 19 in Fig. 1 


P 


3450.13 


(121) 


[414]- 


(120) [515] 


10 549 


p 


47.0 


5.6E-18 




Would be blended with 19 in Fig. 1 


q 


3449.94 


(003) 


[441]- 


(002) [542] 


11468 


p 


38.8 


2.7E-18 


Yes 


Coincides with a weak feature in Fig. 1 


r 


3449.64 


(003) 


[440]- 


(002) [541] 


11468 


o 


38.9 


8.2E-18 


Yes 


Not seen and J—Ka^ 1 



their analysis on Giotto measurements of the ion tempera- 
tures in the coma of Comet Halley (Lammerzahl et al., 1987) 

and assume the electron temperature is approximately equal 
to the ion temperature. The temperature of the electrons 
rises (Faure et al., 2004), and hence the cross-section for 
e~-H20 collisions declines with distance from the nucleus. 
In addition, away from the inner coma, electron density is 
assumed to fall-off as the inverse square of the nuclear dis- 
tance. Hence the greatest effect of electron excitation will 
be within the neutral-neutral coUisionally-dominated inner 
coma, where Te ~200 K. However, Xie and Mumma (1992) 
show that the rate of e~-H2 collisions is still significant 
at distances of 3x10* km from the rmcleus, where To is > 
10 000 K, which represents energies that are able to excite 
the Sv polyad in the H2O molecule. 

The second process also involves electrons. The recom- 
bination of HsO"*" (formed as a result of dissociation of wa- 
ter) gives rise to H2O and H with 25% efficiency, and the 
fragments carry 6.4 eV of excess energy (Jensen et al., 2000). 
This exceeds the bond energy of H2O, but it is very likely 
that a large fraction of the energy is carried by fast H-atoms, 
leaving the H2O in bound states that arc highly vibrationally 
excited. Recent experimental work by Mann et al. (2008) on 
the dissociation of H3O following charge exchange of HsO"*" 
with Caesium revealed that H2O is produced in excited sym- 



metric/asymmetric stretch modes coupled with low bending 
and rotational excitation. The fact that the majority of the 
H3O''' ions undergoing dissociative recombination were ini- 
tially in ground vibrational states with a rotational tem- 
perature of 20-60 K and that the majority of the available 
energy is partitioned to H2O internal energy would appear 
to have direct parallels with vibrational excitation of H2O 
molecules initially in low J, ground vibrational states in the 
coUisionally-dominated region of the cometary inner coma. 

In the last of the three processes, photolysis of H2O 
produces 0(^D) atoms with a quantum yield of 5% under 
quiet Sun conditions (Huebner et al., 1992). 0(^D) is long 
lived ( 110 s), so it can collide with water in the inner coma 
before radiative relaxation. About 1.97 eV ( 16 000 cni^^) 
is available upon collision, more than enough to excite the 
Sv polyad, and + S polyad. Subsequent vibrational cas- 
cade could feed the observed SH lines, but detailed models 
are needed to evaluate this. An important issue is the bal- 
ance between excitation transfer and dissociative reactions. 
Some laboratory work suggests efficient production of two 
OH fragments rather than excitation transfer (Dunlea and 
Ravishankara, 2004). 

Further work is required in order to gain a more com- 
plete understanding of the production of SH lines. This will 
include more high S/N observations in order to correlate 
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the presence of particular SH lines with parameters such 
as cometary activity and nucleocentric distance (i.e., local 
density) and additional laboratory work to understand and 
quantify the physical processes involved. 
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